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À l’occasion de la sortie de son aide-mémoire pour Python 3 aux éditions Ellipses (Le petit Python, 
mai 2017), Richard Gomez nous propose une récréation amusante sur les groupes. Il s’agit d'écrire 
un programme capable de générer tous les groupes finis à isomorphisme près. Cet article a été 
publié dans la revue Mathématice et sur le site Mégamaths. 


Résumé 

La structure d’un groupe est contenue dans sa table de Pythagore. Si on numérote les 
éléments d’un groupe fini d’ordre n, cette table s’identifie naturellement à une matrice carrée 
d'ordre n à coefficients dans {1,...,n}. Il est clair que dans une telle matrice, chaque ligne et 
chaque colonne est une permutation de (1,...,n). Dans le présent article, une telle matrice est 
appelée sudoku. Il est facile d'écrire un programme en Python générant tous les sudokus d’un 
ordre donné, mais attention : un sudoku quelconque ne définit pas forcément un groupe. En 
revanche, si la loi interne définie par un sudoku est associative, alors on a affaire à un groupe. 
Tester une matrice pour savoir si elle est « associative » n’est pas difficile. Nous proposons 
ici un programme capable de trouver tous les sudokus associatifs d’ordre inférieur ou égal à 
6 (au delà, les calculs prennent trop de temps). Notre programme fait ensuite un tri : il ne 
garde qu’un groupe par classe d’isomorphisme. Au final, on se retrouve avec la liste complète 
des groupes d’ordre au plus 6, à isomorphismes près. 


1 Une caractérisation des lois de groupe 
On rappelle les notions de loi interne (loi de composition interne) et de groupe. 


Définition 1. Soit E un ensemble. Une loi de composition interne dans E, est une application de 
Ex E vers E. L'image de (a,b) est noté a X b ou encore ab (notation multiplicative). 


Définition 2. Un groupe est un couple (G, X) où G est un ensemble non vide et X une loi de 
composition interne dans G vérifiant les trois propriétés ci-dessous. 

— Pour tous a,b,ceG@, (ab)c=a(bc). 

—  Îlexiste eE G tel que pour toutaeG,ae=ea—=a. 


Pour tout ae G, il existe a! E G tel que aa'=a'a=e. 


Une loi vérifiant la première propriété est dite associative. L'élément e exhibé dans la deuxième 
propriété s'appelle élément neutre de G. L'élément a’ de la troisième propriété est appelé symé- 
trique de a. On invite le lecteur à voir ou revoir les propriétés fondamentales des groupes dans [1] 
ou [4], par exemple. Nous donnerons dans cette section une condition nécessaire et suffisante pour 
qu’une loi interne soit une loi de groupe. 


Définition 3. Une structure associative est un couple (E, x) où E est un ensemble non vide et 
x une loi de composition associative de E. 


Une structure associative n’est pas un groupe : l’associativité est nécessaire, mais loin d’être 
suffisante. 

Exercice 1. Donner des exemples de structures associatives qui ne sont pas des groupes. 
Définition 4. Une structure sudoku est un couple (E, x) où E est un ensemble non vide et x 
une loi interne dans E vérifiant les deux propriétés suivantes. 

— Pour tous a, bE E, il existe un unique x € E tel que ax =D. 


— Pour tous a,bE E, il existe un unique x E E tel que xa—=b. 


2 SECTION 2 


On notera que dans une structure sudoku, chaque élément est régulier à gauche, autrement dit, 
une égalité a x = a y implique l'égalité x = y (le lecteur le vérifiera par lui-même). De même, chaque 
élément est régulier à droite, ce qui signifie qu’une égalité x a = ya implique x = y. La propriété 
qui suit affirme que la notion de groupe équivaut à la notion de structure sudoku associative. 


Proposition 5. Soit E un ensemble muni d’une loi interne X. Alors (E, X) est un groupe si et 
seulement si (E,x) est une structure sudoku associative. 


Démonstration. Tout groupe est clairement une structure sudoku associative (on invite néan- 
moins le lecteur novice à le vérifier). Il s’agit donc d'établir la réciproque. Soit (E, x) une structure 
sudoku associative. 
1. Soit a€ E. Soitee E tel que 
ae=a 
Un tel élément existe puisque E est une structure sudoku. Montrons que e est neutre à 
gauche. Soit x € E. De ae — a on déduit aex=a x, c’est-à-dire a (e x) = a x, d’où, par 
régularité, ex = x. 
2. Soit f € E tel que fa= a. Montrons que f est neutre à droite. Soit ze E. De fa=a,on 
déduit x fa= x a, c’est-à-dire (x f)a=x a, d’où, par régularité, x f = x. 
3. Montrons que e= f. Puisque f est neutre à droite, e f =e. Puisque e est neutre à gauche, 
ef = f. On a donc e= f. Nous avons prouvé que E possède un élément neutre noté e. 
4. Montrons que tout élément possède un symétrique. Soit a € E. Il existe alors a’ € E tel que 
a a/—=e (nous sommes dans une structure sudoku). Montrons que a” est aussi symétrique à 
gauche de a. De a’e—a” et e= aa, on déduit a'aa'=e a, c’est-à-dire (a’a)a'=e a”, d’où 
a'a=e. 
5. La loi x est supposée associative, et nous avons montré qu’elle vérifie les deuxième et 
troisième axiomes de la définition d’un groupe. Il s’ensuit que (E, x) est un groupe. 


Exercice 2. Soit (E, x) une structure associative. Montrer que (E, x) est un groupe si et seulement si il vérifie 
l’axiome du sandwich : pour tous a,b,ce E, il existe un et un seul x € E tel que aæb=c. Pour en savoir plus 
sur cet axiome, on consultera [7]. 


2 Sudokus associatifs 


Soit E un ensemble fini de cardinal n > 0 muni d’une loi x interne. On numérote les éléments de 
ÆE de sorte que ce dernier s’écrive 


E={a1,.….,an} 
On peut maintenant associer de manière naturelle une matrice P d’ordre n au couple (E, x). En 
effet, soit (i,j) € {1,...,n}?. On note k l’élément de {1,...,n} tel que 
di a j — x 
On pose alors 
P;5=k 


La matrice P ainsi définie est à coefficients dans {1,...,n}. Bien entendu, P dépend de la manière 
dont on a numéroté les éléments de ÆE. Néanmoins, une fois numérotés les éléments de E, il y a 
correspondance bi-univoque entre les lois internes de E, et les matrices d’ordre n à coefficients dans 


{Lis 


Supposons maintenant que (E, X) est une structure sudoku avec éléments numérotés. On note P 
la matrice associée. Le fait que chaque équation a x = b d’inconnue x possède une et seule solution 
dans FE, implique que chaque ligne de P est une permutation de (1,...,n). De même, le fait que 
chaque équation x a = b possède exactement une solution, implique que chaque colonne de P est 
une permutation de (1,...,n). Dans une telle matrice, aucune ligne ne se répète (aucune colonne 
non plus). Nous allons nous intéresser à ce type de matrices (pour simplifier notre programme 
Python à venir, nous ne nous intéresserons qu’au cas où la première colonne et la première ligne 
correspondent à la permutation triviale). 
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Définition 6. Soit P une matrice carrée d'ordre n à coefficients dans {1,...,n}. On dit que P est 
un sudoku si chaque ligne et chaque colonne est une permutation de (1,...,n), la première ligne et 
première colonne étant égales à (1 … n). 


On aura aussi besoin de la définition ci-dessous : 


Définition 7. Soit P une matrice d’ordre n à coefficients dans {1, …, n}. On dit que P est 
associative si la loi définie par P est associative. 


La proposition ci-dessous est une simple conséquence de la proposition 5 de la section précé- 
dente : 


Proposition 8. Soit G un groupe fini. On numérote les éléments de G de sorte que le neutre porte 
le numéro 1. Alors la matrice associée à cette numérotation est un sudoku associatif. Réciproque- 
ment, si P est un sudoku associatif d'ordre n, alors il définit une loi de groupe sur {1,...,n} pour 
laquelle 1 est l'élément neutre (tout sudoku associatif vient d’un groupe). 


Ainsi, pour trouver tous les groupes d’ordre n (à isomorphisme près), il suffit de construire tous 
les sudokus d'ordre n, et ne garder que ceux qui sont associatifs. Sans oublier de ne garder qu’un 
exemplaire par classe d’isomorphisme. C’est ce que nous ferons avec Python. 


Remarque 9. D'un point de vue strictement mathématique, la correspondance entre groupes 
finis et sudokus associatifs n’a pas grand intérêt. Chercher à classifier les groupes finis en cher- 
chant les sudokus est une démarche grossière. De plus, cette méthode devient inefficace dès que 
l’ordre devient grand. Pour une recherche intelligente (et intéressante), on consultera [1] et [3], par 
exemple. Le but du présent article est modeste : écrire un programme simple en langage Python 
capable de trouver tous les groupes d’ordre au plus 6 (à isomorphisme près) en un temps très court. 


Exercice 3. Vérifier qu’il n’existe qu’un sudokus d’ordre 3. Pourquoi ce sudoku est-il forcément associatif ? 
Combien y a t-il de sudokus d’ordre 4 ? 


Exercice 4. Expliquer pourquoi le fait d’avoir restreint la définition des sudokus aux matrices dont la première 
ligne et première colonnes coïncident avec (1 … n) (définition 6) rend la proposition 5 inutile pour démontrer 
la proposition 8. 


3 Python et permutations 


Pour détecter les éventuels isomorphismes entre groupes, notre programme aura besoin de générer 
toutes les permutations d’un degré donné. Dans cette section, nous montrons plusieurs manières 
d'y parvenir. Bien sûr, la liste des méthodes proposées ici n’est pas exhaustive. Le lecteur pressé 
peut sauter cette section. 


3.1 Convention 


On rappelle qu’une permutation d’un ensemble E est une bijection de Æ sur lui-même. L'ensemble 
des permutations de Æ muni de la loi de composition des applications est un groupe noté S(E). 
Le groupe des permutations de {1,...,n} est appelé groupe symétrique de degré n, et se note Sa. 
Donnons un exemple. La permutation de {1,2,3} définie par 

1 3 

21 1 


3 — 2 


1273 
3 1 2 


(les antécédents en haut et les images en bas). En informatique, on encode cette bijection avec 
Puplet (3,1,2), tout simplement. Nous ferons de même. 


se note 
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Exercice 5. Lister toutes les permutations de {1,2,3} (en suivant notre convention). 


3.2 Algorithme par itérations 


Nous proposons ici un programme bestial (et peu efficace) pour calculer l’ensemble des permuta- 
tions de {1,2,3}. Il consiste à générer tous les éléments de {1,2,3}% et ne conserver que les uplets 
sans répétitions. 


On voit mal comment convertir ce code en une fonction capable de calculer l’ensemble des per- 
mutations d’un itérable de longueur inconnue a priori. Néanmoins, l’idée exposée ici peut nous 
dépanner dans des situations particulières (quand il s’agit d’un besoin très simple). 


Remarque 10. Les objets de type set, comme par exemple {1,2,3}, sont étudiés au chapitre 23 
de [5]. Les objets de type list, comme PERM, au chapitre 20. Les objets de type tuple, comme 
(x,y,z), au chapitre 21. La fonction len renvoie la longueur de n'importe quel objet itérable, 
de sorte que len({x,y,z}) renvoie le cardinal de l’ensemble {x,y,z}. La ligne PERM = PERM + 
[(x,y,z)] peut être remplacée par PERM.append((x,y,z)). 


3.3 Fonction itertools.permutations 


Si on envoie un itérable L à la fonction permutations du module itertools, cette dernière 
retourne un iterator générant toutes les permutations de L. Importons permutations depuis le 
module itertools : 


Créons un itérable nommé L (une liste par exemple) : 


Pour avoir accès aux éléments de permutations(L), on peut par exemple itérer avec une boucle 
for : 


On retrouve bien les 6 permutations de [1,2,3]. Chacune d'elles est encodée par un tuple. 


Remarque 11. Les modules (et les imports) sont décrits au chapitre 5 de [5]. Les listes au 
chapitre 20. Les uplets (tuple) au chapitre 21. Les itérators (generators) et le module itertools 
au chapitre 22. 
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On peut donner à permutations une liste contenant des répétitions : 


Dans tous ces exemples, les items de l’itérable envoyé à permutations étaient des nombres entiers, 
mais ceci n’est nullement une obligation. On peut envoyer une liste de chaînes de caractères, ou 
une liste d’items de n’importe quelle nature. 


3.4 Une fonction récursive 


Dans cette section, nous écrivons nous-mêmes une fonction calculant les permutations. Son utilisa- 
tion est simple : on lui donne un itérable, et elle retourne la liste des permutations de cet itérable. 
Chaque permutation est encodée par une liste. 


La première instruction équivaut à ceci : si L est vide, on retourne [] (liste vide). La deuxième 
instruction équivaut à ceci : si L ne possède qu’un item, on retourne [list(L)], liste contenant 
l'unique permutation de L (la permutation triviale, encodée par une liste). La troisième instruction 
est récursive (voir plus bas). 


Remarque 12. La structure if...elif...else est décrite au chapitre 9 de [5]. Les booléens 
sont décrits au chapitre 17 (on comprendra alors l'instruction if not L). La syntaxe pour écrire 
une fonction (def, return, etc.) est décrite au chapitre 13. 


La récursivité écrite ici est basée sur une idée simple. Regardons un exemple : pour écrire toutes 
les permutations de {1, 2, 3}, on commence par les permutations (1, x, *), ensuite on passe aux 
(2, x, x), et on termine avec les (3, *, x). Pour obtenir tous les (1, x, x), on concatène (1) avec 
chaque permutation de {2,3}. Pour obtenir les (2,*,*), même chose : on concatène (2) avec chaque 
permutation de {1, 3} et ainsi de suite. Tout le monde voit la relation de récursivité : chaque 
permutation de (1,2,3) s'écrit en concaténant chaque (a) avec chaque permutation de {b,c}, ce 
dernier étant l’ensemble obtenu en retirant a de {1, 2, 3}. Cette concaténation s'écrit [L[i]]+P 
dans notre script, où L[i] décrit tous les items de L, et P décrit toutes les permutations de la liste 
«L privée de Li] ». Pour que Li] décrive tous les items de L, il suffit que i décrive toutes les 
positions possibles dans L, c’est ce que nous provoquons en écrivant for i in range(len(L)). 
Pour que P décrive toutes les permutations de la liste L privée de L[i], nous avons écrit for P in 
perm(Ll:i]+L{i+1:]). Le lecteur vérifiera que L[:i]+L[i+1:] est bien la liste L privée de son 
item [il]. 


Remarque 13. Le type range est décrit au chapitre 16 de [5]. La fonction 1en retourne la longueur 
d’un itérable. Les syntaxes L[:i] et L[i:] sont décrites au chapitre 20 (consacré aux listes). 


On observera que notre fonction retourne bien la liste de tous les LTi]+P attendus grâce à une 
écriture typique de Python : l'écriture en compréhension. Schématiquement, la fonction retourne : 


[ [L[i]]+P for i in here for P in there |] 
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Remarque 14. L'écriture d’une liste en compréhension est décrite en détail aux sections 20.9 et 
20.10 de [5] (chapitre 20 consacré aux listes). 


La récursivité s’arrête forcément dans cette fonction. En effet, l'appel perm(L) déclenche des appels 
perm(Ll:i]+Lli+1:]) où la liste L[:i]+L{[i+1:] est strictement plus courte que la liste initiale. 
On finit donc par tomber sur des listes de longueur 1, c’est-à-dire le cas traité dès la deuxième 
instruction (arrêt de la récursivité). 


Exercice 6. Exécuter à la main la fonction perm pour L=[1,2,3]. 


3.5 Un algorithme de Donald Knuth 


On trouve dans [2] un algorithme basé sur la notion d'ordre lexicographique. L'idée consiste à 
munir chaque groupe $, de cette relation d'ordre. En vertu de la convention faite à la section 3.1, 
l’ensemble S3, par exemple, s'écrit 


S3= {(1,2,3),(1,3,2),(2,1,3),(2,3,1),(3,1,2),(3,2,1)} 


et comparer lexicographiquement (2,1,3) avec (3,1,2) revient à comparer les entiers 213 et 312. On 
trouve bien sûr (2,1,3) <(3,1,2). Le lecteur attentif aura remarqué que nous avons listé les éléments 
de S3 dans l’ordre croissant. Le noyau dur de l'algorithme de Knuth est une fonction capable de cal- 
culer la permutation qui arrive juste après une permutation donnée, selon l’ordre lexicographique. 
Cette fonction, nous l’appellerons suivant, de sorte que suivant ([2,1,3]) retournera [2,3,1]. 
L’argument de cette fonction est obligatoirement une liste (ie. une donnée de type list). Voici la 
traduction en langage Python de la fonction suivant : 


Remarque 15. La fonction suivant est intéressante du point de vue mathématique et algorith- 
mique. Le lecteur pourra l’étudier dans [2]. Son code Python en revanche ne présente pas grand 
intérêt, si ce n’est la nécessité dès la première ligne de faire une copie de la liste L avec la méthode 
copy, afin d'éviter un effet de bord. Ces derniers sont décrits à la section 13.13 du chapitre 13 de 
[5]. La méthode copy et le phénomène d’aliasing sont largement documentés dans [5]. 


Testons la fonction suivant : 
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On peut même demander quelle est la « permutation » de (1,1,1,3) — on nous pardonnera cet 
abus de langage — qui vient juste après (1,1,3, 1) : 


Analysons ce que fait la fonction suivant pour construire suivant ([10,3,4,9,1,5,8,7,6,2]). 
D'abord, elle cherche deux items consécutifs a ,b tels que a<b, en commençant par la fin 6,2]. Elle 
trouve 5,8 : 


10[3[4]9/1]5]8[7[6[2 


Ensuite elle cherche l’item strictement supérieur à notre 5 grisé sur la figure, en commençant par 
la fin, 2. Elle trouve 6 : 


10]3[4]911]5[8]7[6[2 


En position [i] nous avons 5, et en [j] nous avons 6. La fonction permute ces deux items, 


10]3[4[911]6[8[7[5[2 


augmente [i] d’un cran, positionne [j] à la fin : 


10[3[4]911[6[8[7[5]2 


et entame une série de permutations selon le principe suivant : permutation des deux items gris, 
augmentation de Li], diminution de [j], permutation des deux gris, et ainsi de suite jusqu’à ce 
que les deux items gris se rejoignent (ou se croisent) : 


10]3[4[9]1]6]2]7]5]8 


101314191116/21715|8 


10]3[4[9]1]6[2]5[7[8 


La fonction retourne (10, 3,4,9,1,6,2,5,7,8). 


Maintenant que nous avons compris à quoi sert la fonction suivant, nous pouvons écrire un 
algorithme utilisant cette dernière pour construire la liste des permutations d’une liste donnée : 


Si les items de L sont dans l’ordre croissant, l'appel knuth(L) retourne la liste des permutations 
de L : 


Si les items de L ne sont pas dans l’ordre croissant, on se doute du résultat : 
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Exercice 7. Si L contient la liste vide [], les appels perm(L) (voir section 3.4) et knuth(L) ne retournent pas 
la même chose. D'un point de vue mathématique, laquelle a raison ? 


Exercice 8. Réécrire la fonction knuth afin que knuth(L) retourne la liste des permutations de L même si cette 
dernière n’est pas donnée dans l’ordre croissant. Indication : utiliser la fonction sorted ou la méthode sort. 


On peut modifier knuth pour en faire un generator, à l'instar de itertools.permutations : 


Dans ce cas, knuth retourne non pas une liste mais un objet assez complexe (un generator iterator) 
qui a vocation à être itéré avec une boucle par exemple : 


Remarque 16. La dernière version de knuth exploite des spécificités fortes de Python telles 
que yield et raise. L’instruction yield est étudiée au chapitre 22 de [5] (chapitre consacré aux 
générateurs), et l'instruction raise au chapitre 12 (gestion des exceptions avec try-except-else). 


3.6 Algorithme par insertions 


Nous nous contenterons de décrire l’algorithme par insertions sur un exemple simple : le calcul 
de l’ensemble des permutations de {1,2,3,4}. Pour ce faire, l'algorithme commence par écrire les 
permutations de {1}, puis celles de {1, 2}, puis celles de {1, 2, 3}, pour enfin produire celles de 
{1,2,3,4). 


e  Permutations de {1} : (1). 


e  Permutations de {1,2} : on insère 2 dans chaque permutation de {1}. On trouve : (2,1), 
(1,2). 
e  Permutations de {1,2,3} : on insère 3 dans chaque permutation de {1,2} : 
— (2,1) donne (3,2,1), (2,3,1), (2,1, 3). 
— (1,2) donne (3,1,2), (1,3,2), (1,2, 3). 
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e  Permutations de {1,2,3,4} : on insère À dans chaque permutation de {1,2,3}. 
— _(3,2,1) donne (4,3,2,1), (3,4,2,1), (3,2,4,1), (3,2,1,4). 
— etc. 


Exercice 9. Traduire l’algorithme par insertions en langage Python. 


4 Programme pour une recherche exhaustive des groupes 


Dans cette section, une matrice d'ordre n est une matrice carrée n x n à coefficients dans l’ensemble 
{0,...,n—1}. Nous aurions préféré travailler avec des coefficients dans {1,...,n}, mais nous avons 
préféré nous adapter à Python qui, comme chacun sait, numérote les items d’une liste en partant de 
0 et non de 1. Nous encodons une matrice d'ordre n avec une liste de listes. Par exemple, la matrice 


112 
M=| 110 

2 1 2 
est encodée par : 


Le coefficient grisé s'appelle normalement M3 (ligne 2, colonne 3), mais nous préférons nous 
adapter à Python en l’appelant M2 : 


Nous regarderons M comme la table de Pythagore d’une loi interne de {0, 1, 2} notée x. Par 
exemple, ici nous avons : 


1xX2—=0 


4.1 Fonction générant tous les sudokus d’ordre n 


La fonction sudoku présentée ici reçoit un entier naturel non nul n en argument, et retourne la liste 
de tous les sudokus d’ordre n. Montrons sur l’exemple n=4 comment cette fonction travaille. En 
fait, elle fabrique un arbre généalogique. La première génération est (0 1 2 3 ), ancêtre commun de 
tous les sudokus d’ordre 4. Pour construire les enfants de ce dernier, on calcule les lignes acceptables 
(celles qui conservent la propriété d’être un sudoku, nous montrerons plus loin comment). On 


en trouve trois : (1 0 3 2), (1 2 3 0) et (1 3 0 2). On peut donc dire que (0 1 2 3) enfante 
0 1 2 3 0 1 2 3 ü. 1 23 : 

( rs }: ( mt ) et ( à di D ). On recommence avec chacun de ces trois presque-sudokus. 

Voici le résultat final : 


(0 123) 


0123 0123 0123 

1032 1230 1302 
0123 0123 123 0123 
1032 1032 2 3 0 1302 
2301 2 3 1 0 3 01 2031 
0123 0123 0123 0123 
1032 1032 1230 1302 
2301 2.31 0 2301 2031 
3 21 0 3 201 3 012 3 2 1 0 
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Nous voyons que l’appel sudoku (4) retournera quatre sudokus d'ordre 4. Montrons maintenant 
comment on calcule les lignes qu’un presque-sudoku peut accepter. Prenons par exemple 


0123 4 
a=( 13042 ) 
Il s’agit bien d’un presque-sudoku. Ici aussi, la méthode consiste à fabriquer un arbre généalogique. 
La première génération est la presque-ligne (2). Passons à la deuxième génération. Après le 2, on 
ne peut pas mettre un 2. De plus, la deuxième colonne de À contient déjà 1 et 3, donc les seuls 
coefficients acceptables sont 0 et 4. Autrement dit, les enfants de (2) sont (2 o) et (2 4); c’est 
la deuxième génération. On passe à la troisième génération. D’après le même raisonnement, (2 0) 


donne (2 0 1)},(2 0 3)et(2 0 4), tandis que (2 4) donne (2 4 1)et(2 4 3). On continue ainsi 
de suite jusqu’à la cinquième génération : 


D — 


(2 0) (2 4) 


a 


(201) (203) (204) (241) ( ) 


| | | _…—.._ 


(SSL (2040) CAT) (fat (AXEL 


(ÉDELTSIN Ed Ts) Cid Ir ds To) 


ND 
FIN 
[ee] 


Ainsi, les lignes que À peut accepter sont (2 0 4 1 3),(2 4 1 0 3),(2 4 3 0 1)et(2 4 3 1 0). 
On en déduit les enfants de À : 


Voici le code de notre fonction sudoku : 


PROGRAMME POUR UNE RECHERCHE EXHAUSTIVE DES GROUPES IT 


Variables utilisées : 
— GC: génération courante (contient les presque-sudokus de la génération en cours). 


— _GN: génération nouvelle. C’est la génération venant juste après GC. Chaque presque-sudoku 
M de GC enfante une liste DM (c’est la liste des presque-sudokus dérivant de M). Aïnsi, pour 
chaque M de GC, on enrichit GN avec les éléments de DM grâce à l'instruction GN = GN + DM. 


— _gc: génération courante pour les lignes (les presque-lignes de la génération en cours). 


— _gn: génération nouvelle. C’est la génération venant juste après ge. Chaque presque ligne 
L de gc enfante une liste DL (c’est la liste des presque-lignes dérivant de L). Ici aussi nous 
retrouvons naturellement l’instruction gn = gn + DL. 


Remarque 17. Les listes et la fonction list font l’objet du chapitre 20 de [5]. Les boucles for font 
l’objet du chapitre 10. Les ensembles et la fonction set font l’objet du chapitre 23. La définition 
d’un ensemble en compréhension comme par exemple {M[i][j] for i in range(0,i)} y est 
expliquée en détail. 


4.2 Alternative bestiale pour sudoku 


Au lieu de calculer intelligemment l’ensemble des lignes acceptables pour un presque-sudoku donné, 
on peut tout simplement générer toutes les lignes possibles, c’est-à-dire toutes les permutations de 
(0 1 … n—1), et les tester une par une pour savoir lesquelles on garde. Soit M un presque-sudoku 
de dimensions m x n (m<n). Soit L une ligne de longueur n (c’est-à-dire une matrice 1 x n). On 
écrit ici une fonction retournant True si l’ajout de L à M conserve la propriété d’être un sudoku, 
False sinon : 


M 


L 


Il est clair que ZL est acceptable si et seulement si pour tout 7, le coefficient L; ne se trouve pas 
dans la colonne j de M. Voici la fonction : 


La fonction bestiale remplaçant sudoku serait alors : 
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On notera que sudoku_bestial utilise la fonction permutations du module itertools. 


Remarque 18. Les méthodes ordinaires appartenant au type list, comme par exemple la 
méthode remove appliquée à la liste nommée support dans ce code (ligne 6), sont décrites à 
la section 20.13 de [5]. 


4.3 Fonction pour savoir si une matrice est associative 
Dire qu’une matrice M d’ordre n est associative, signifie que la loi x qu’elle définit sur {0,..,n—1} 
vérifie la propriété 
GDE=I GR) 
autrement dit que 
Mu,;k= Mi, 


is 


Voici le code de notre fonction testant l’associativité : 


Une fois construits tous les sudokus d'ordre n, nous les testerons un par un pour ne garder que les 
associatifs. On rappelle que chaque sudoku associatif représente un groupe. 


Remarque 19. Les opérateurs de comparaison comme != (qui signifie Æ) sont étudiés à la section 
17.2 du chapitre 17 (booléens) de [5]. 


4.4 Fonction pour savoir si une application est un isomorphisme de 
groupes 


Soient M et N des sudokus associatifs d'ordre n. Chacune de ces matrices est considérée comme la 
table de Pythagore d’une loi de groupe sur {0,...,n—1}. On note M et NW les groupes définis par M 
et N respectivement. Soit f une bijection de {0,...,n—1}. On écrit ici une fonction isomorphisme 
répondant True si f: M —N est un isomorphisme de groupes, et False sinon. Il s’agit de vérifier 
que l’on a 


f(ab)= f(a) f(b) 
pour tous a, bE{0,..,n—1}. Autrement dit, vérifier que l’on a 
FMav) = N j(a).f(b) 


Dans notre programme, la fonction isomorphisme ne traitera que des bijections f vérifiant f(0)—0. 
On rappelle que les premières ligne et colonne d’un sudoku sont par définition triviales, de sorte que 
tout sudoku associatif est la table d'une loi de groupe sur {0,...,n — 1} où 0 joue le rôle d’élément 
neutre. Il est donc inutile de tester l'égalité f(ab) = f(a) f(b) lorsque a ou b est nul. 
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Remarque 20. D'un point de vue mathématique, £ est une bijection, mais d’un point de vue 
informatique, f est un tuple (nous rappelons que chaque permutation est encodée ici avec un 
uplet), c’est la raison pour laquelle dans notre script, une image f(x) s'écrit f [x], et non pas f (x). 


4.5 Fonction pour savoir si deux groupes sont isomorphes 


Nous reprenons les notations de la sous-section précédente. Nous proposons un algorithme bestial 
pour savoir si M et NW sont isomorphes. La fonction isomorphe regarde s’il existe une bijection 
f:M = N fixant 0 qui soit une morphisme de groupe. 


La fonction permutations appelée à la troisième ligne est celle qui est offerte par le module 
itertools (section 3.3). Il faudra penser à l’importer. 


Remarque 21. Les objets de type tuple, comme le f et le (0,) de ce code sont étudiés au chapitre 
21 de [5]. D'un point de vue mathématique, un objet de type tuple est un uplet. L’instruction 
return f (ligne 6) aura le même effet qu’un return True. Pour comprendre pourquoi, on pourra 
consulter le chapitre 17 de [5] consacré aux booléens. 


4.6 Fonction pour savoir si un groupe donné est déjà représenté dans la 
liste 


Chaque fois que nous générons un groupe, nous regardons si notre stock ne contient pas déjà un 
groupe qui lui serait isomorphe. La fonction est_dans reçoit deux arguments : un sudoku associatif 
M et une liste de sudokus associatifs L. 


4.7 Fonction générant tous les groupes à isomorphisme près 


La fonction groupes reçoit un entier naturel n non nul en argument, et retourne la liste de tous 
les groupes d'ordre n à isomorphisme près. Ainsi, chaque classe d’isomorphisme possède un repré- 
sentant et un seul dans la liste retournée. 
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Remarque 22. Quand Python calcule le booléen associative(M) and not est_dans(M,G) lors 
de l'instruction if de la ligne 5, il commence par calculer associative(M). Si celui-ci vaut False, 
il arrête le calcul car il sait que le booléen cherché vaut False lui aussi. Si tel n'avait été le cas, 
nous aurions décomposé ce simple if en un double if afin de gagner en vitesse d'exécution. Nous 
décrivons la manière dont Python utilise les opérateurs logiques tels que and et or aux sections 
17.5 et 17.6 de [5]. 


5 Résultats de notre recherche 


Nous ne parlerons pas du cas n = 1. Nous avons utilisé le décorateur décrit à la section 14.7.3 du 
livre [5] pour chronométrer notre programme. 


5.1 Groupes d’ordre 2 
— Recherche des sudokus : 0,000 02 secondes (1 seul sudoku). 
— Recherche des groupes : 0,000 03 secondes (1 seul groupe). 


Groupe trouvé : ( : À }: On reconnaît Z/27Z. C'est conforme à la théorie, Cf [6]. 


Remarque 23. Nous avons chronométré le temps utilisé pour trouver les sudokus, puis le temps 
utilisé pour trouver les groupes parmi ces derniers. 


5.2 Groupes d’ordre 3 
— Recherche des sudokus : 0,000 05 secondes (1 seul sudoku). 
— Recherche des groupes : 0,000 03 secondes (1 seul groupe). 
1 2 
Groupe trouvé : | 2 0 | 
201 


On reconnaît le groupe Z/3Z. C’est conforme à la théorie. 


5.3 Groupes d’ordre 4 
— Recherche des sudokus : 0,000 1 secondes (4 sudokus, tous associatifs). 


— Recherche des groupes : 0,000 1 secondes (2 groupes). 


0 1 2 3 0 1 2 3 
su | À 03.2 1032 

Groupes trouvés : Lana lee so 
3 2 10 3201 


On reconnaît les groupes (Z/2 Z)? et Z/4 7%, respectivement. C’est conforme à la théorie. Voici 
une activité permettant de comprendre comment on a fait pour reconnaître Z /4 Z. On note 
Z/4Z={0,1,2,3}. La table de ce groupe est alors 


On compare les deux tables avec la fonction isomorphe : 
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Python nous répond en donnant un isomorphisme. 


5.4 Groupes d’ordre 5 
— Recherche des sudokus : 0,002 secondes (56 sudokus, parmi eux 6 associatifs). 


— Recherche des groupes : 0,000 7 secondes (1 groupe). 


0 1 2 3 4 
12340 
Groupe trouvé : | 23 401 
3 4 0 1 2 
4 0 1 2 3 


On reconnaît le groupe Z /5 Z. C’est conforme à la théorie. Voici une activité permettant de 
comprendre comment on à fait pour savoir que parmi les 56 sudokus trouvés, il y en à 6 associatifs : 


Note. Le lecteur a compris que les 6 sudokus associatifs trouvés représentent le «même» groupe. 


5.5 Groupes d’ordre 6 


— Recherche des sudokus : 0,75 secondes (9408 sudokus, parmi eux 80 associatifs). 


TROT 


— Recherche des groupes : 0,12 secondes (2 groupes). 


Groupes trouvés : et 


H © ND H © 
& Où ND © © H 
k © à X ND 

kH & Où ND 

DhRhOoOumR 

D © H À ot 
a & © DR © 
WDABOH 
& hr © © ND 
OR ND & & 


5 0 3 2 
Le lecteur vérifiera que le premier est Z/6Ze 


+ Hop un w 
D Doors 


le deuxième est S3 (groupe symétrique de degré 3). 
5.6 Groupes d’ordre 7 


Nous n'avons pas eu la patience d’attendre la fin de l'exécution. La fonction sudoku nécessite 9 
minutes pour construire la quatrième génération de presque-sudokus. Les chiffres sont spectacu- 
laires : 


— deuxième génération : 309 presque-sudokus ; 
— troisième génération : 35 792 presque-sudokus ; 


— quatrième génération : 1 293 216 presque-sudokus. 


6 Programme pour une pêche aléatoire 


Pour une pêche aux groupes basée sur la chance, nous aurons besoin d’une fonction générant 
aléatoirement une permutation : 
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Attention : cette fonction utilise la fonction shuffle du module random. On n'oubliera pas de 
l’importer au début du script : 


Remarque 24. Le module random est étudié en détail à la section 7.11 du chapitre 7 (calcul et 
arithmétique) du livre [5]. 


Nous aurons ensuite besoin d’une fonction générant aléatoirement un sudoku : 


Nous aurons enfin besoin d’une fonction générant successivement des sudokus jusqu’à tomber (par 
hasard) sur un sudoku associatif (ie un groupe). Cette fonction se comporte comme un pêcheur 
sélectif. 


Remarque 25. L’instruction while 1 équivaut à while True (boucle infinie). Le chapitre 17 
de [5] explique comment Python transforme automatiquement n'importe quel type de donnée (en 
l'occurence ici l’entier 1) en un booléen. 


Nous proposons d'organiser la pêche comme ceci : 


Quelques explications : à la fin de l’exécution, S[n] contient tous les groupes d'ordre n pêchés 
(maximum un groupe par classe d’isomorphisme). Le programme ne traite que les cas n € {2,3,4, 
5,6}. Pour chaque n, on appelle 30 fois la fonction groupe_alea. À chaque appel, on regarde si 
le groupe obtenu est nouveau ou pas. S'il est nouveau, on le stocke, sinon, on l’abandonne. 


Voici le produit de notre pêche (juste une exécution) : 


—  S[2] a retenu ( . : } c’est-à-dire Z/27: 
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1 2 
—  S[3] a retenu | 1 2 0 |, c’est-à-dire Z/37Z; 
20 1 
ds :235 0 1 23 
_ 1 302 1032 ja RL 2. 
S[4] a retenu montée , c’est-à-dire Z/47Z et (Z/27Z)°: 
3210 3 2 1 0 
0 1234 
12340 . 
—  S[5] à retenu | 2 3 4 0 1 |, c’est-à-dire Z/57; 
3 4 0 1 2 
40123 
012345 012345 
134250 105432 
. 240513 2 3 0 1 5 4 act A A 
S[6] a retenu na ES Poe: , c'est-à-dire Z/6Z et Sa. 
4510 3 2 453201 
5 0 3 1 2 4 541023 


Le hasard a voulu que nous attrapions tout ce qui existe! Voici les durées enregistrées lors de cette 
exécution (arrondies au millionnième de seconde) : 


durée en secondes 
0,000 332 
0,000 942 
0,004 602 
0,102 600 
6,554 921 


SOURIS 


Il est intéressant de noter combien de fois la fonction groupe_alea doit appeler sudoku_alea pour 
tomber sur une loi associative (un groupe ici). Une petite modification de notre programme nous 
a permis d'obtenir les chiffres ci-dessous (sur une exécution) : 


Nombre d’appels en moyenne 
1 
1 
1 
&6, 6 
&115,7 


ol ou &|w| ns 


On constate que pour obtenir un sudoku associatif d'ordre 6, il faut, en moyenne, générer 
116 sudokus. On rappelle que pour une exécution de notre programme de pêche, l’appel 
groupe_alea(6) a lieu 30 fois, c’est donc une moyenne sur 30 mesures. 


7 Perspectives 


Il y a énormément de choses à améliorer dans nos scripts. On peut les rendre plus rapides afin 
de pouvoir espérer traiter quelques cas au-delà du cas n —6 (quelques jours après avoir rédigé cet 
article, l’auteur a résolu le cas n—7). On notera par ailleurs que notre fonction sudoku_alea est 
incapable de trouver un sudoku d’ordre 13 en un temps raisonnable. Nous ne savons même pas 
si elle en est capable tout court ! Une activité intéressante, par exemple, consisterait à écrire une 
fonction permettant à coup sûr de générer un sudoku d'ordre 13 (ou plus) en s'inspirant d’un 
algorithme de type backtracking search, par exemple. Il y a beaucoup à faire... 
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